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ABSTRACT 


Let  g^  (i  *  1,  . ..,  m)  be  smooth  vector  fields  on  Rd,  and  let  Tn_1  , 

their  Taylor  expansions  of  order  n  -  1  at  the  origin.  The  system 
m 

x(t)  *  J  g^xft)  )u^(t) ,  x(o)  -  0  e  R^  generates  an  input-output  map 


i-1 

u( • )  ♦  x(u,* )  whose  n-th  order  Taylor  approximation  xr(u,»)  can  be 

obtained  by  computing  the  n-th  Picard  iterate  for  the  reduced  system 

x(t)  *  (T^HxIt))  u^t),  x(o)  *  0,  discarding  the  terms  of  order 
i-1 

>  n.  For  z  e  R?,  directional  error  bounds  of  the  form 


<  z,  x(Uft)  -  xn(u,t)  >  |  <  C  t”~p  J*  |  J®  gi(o)ui(o)  dc|pds 

°  i-1 


can  be  given.  These  estimates  improve  those  supplied  by  the  classical 
Taylor's  theorem  and  yield  results  concerning  local  non-controllability. 
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SIGNIFICANCE  AND  EXPLANATION 


Consider  the  control  system  on 


R 


A, 


x(t )  ■  g. (x(t) )*u. <t) ,  x(o)  *  0. 
*  i-1  1  1 


Given  a  control  u«  in  general  there  exists  no  explicit  formula  to  exactly 

compute  the  corresponding  trajectory  t  ♦  x(u,t).  -The- present-  paper  provides 

a  simple  method  of  constructing  approximations,  xn(u,t)  of  x(u,t),  which  are 

linear  combinations  of  integrals  of  the  controls,  u^  and  which  can  be 

explicitly  computed.  The  coefficients  of  these  linear  combinations  are 

obtained  from  the  partial  derivatives  of  the  vector  fields  g±  at  the  origin, 

using  a  Picard  iterative  procedure . 

Here  the  approximation  error  is  the  vector  eR  (u,t)  ■  xn(u,t)  - 

x(u,t).  New  and  precise  bound#  on  the  length  of  en<u,t),  as  well  as  on  the 

else  of  the  projections  of  e  (u,t)  on  some  special  subspaces  of  fA  are 

n 

given.  „ 


The  present  formula  can  be  used  for  a  systematic  study  of  a  control 
system  over  a  short  Interval  of  time,  in  much  the  same  way  as  the  classical 
Taylor  expansion  is  used  to  study  the  local  properties  of  differentiable 


functions • 

r\ 


\ 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


ON  THE  TAYLOR  APPROXIMATION  OF  CONTROL  SYSTEMS 
Albarto  Bressan 

Conaidar  an  autonoaxma  nonlinaar  ayatea  on  R^  with  control  antaring  linaarlyi 
x(t)  -  f(x(t))  +  I*  g. (x(t) )  u  (t) 

i-1  1  1 

(1.1) 

x(o)  -  0,  t  >  0. 

Systems  of  the  fora  (1.1)  were  firat  studied  in  a  paper  by  Hermes  and  Haynes  [4]  and 
received  continued  attention  since  than.  If  f  and  the  g^ *s  are  smooth  vector  fields, 
then  (1.1)  yields  a  smooth  input-output  map  4»  (.'  ( [0,»)»  iP)  *  C([0,»)|  Rd)  defined  by 

4(u)(t)  -  x ( u . t ) .  where  x(u,«)  is  the  unique  solution  of  (1.1)  corresponding  to  the 
control  u.  In  general  there  exist  no  explicit  formulas  giving  the  trajectory  x(u,») 
directly  in  terms  of  the  control.  It  is  therefore  natural  to  look  for  aomm  computable 
approximation  of  the  map  4.  The  Taylor  expansion  of  4  in  terms  of  Volterra  kernels 
was  considered  by  Brockett  I 1 ] .  The  local  aproximation  of  a  control  system  of  the  form 

i(t)  -  y  g  (x(t) )u. (t) ,  x(o)  -  0  (1.2) 

i-1  1  1 

by  means  of  an  auxiliary  (nilpotent)  system  is  studied  in  (5).  For  analytic  systems, 
expansions  in  formal  power  series  are  given  in  [3]. 

In  the  present  paper  we  approximate  the  input-output  may  4  generated  by  (1.2)  with 
linear  combinations  of  certain  iterated  integrals,  here  called  integral  monomials.  Using 
functional  analytic  techniques,  we  derive  a  simple  procedure  to  recursively  compute  the 
coefficients  of  the  Taylor  expansion  for  4  in  terms  of  the  Taylor  coefficients  of  the 
g^'s  at  the  origin.  No  analyticity  assumptions  are  needed.  Our  main  concern  is  a 
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precise  estimate  of  the  approximation  errors.  Since  the  system  (1.2)  is  linear  in  u, 
these  errors  are  essentially  due  to  the  nonlinearity  of  the  vector  fields  gi#  which 
depend  only  on  x.  We  thus  obtain  bounds  which  depend  not  on  the  size  of  the  control 
u ( • )  but  only  on  the  norm  of  a  first  order  approximation  to  the  trajectory  x(u,<). 

The  estimates  contained  in  Theorems  1  and  2,  in  sections  5  and  6  respectively,  improve 
those  supplied  by  the  classical  Taylor's  theorem  and  are  sharp  enough  to  yield  results  on 
xocal  non-controllability  for  the  system  (1.1).  Examples  are  provided  in  sections  4  to 
6.  Some  related  results  recently  appeared  in  [6] . 

$  2.  Preliminary  Abstract  Results 

1. 

Given  two  Banach  spaces  E  and  F,  k  >  0,  we  denote  by  L  (E,F)  the  space  of 

continuous  k-llnear  mappings  A  from®  E  *  ExEx.-.xE  (k  times)  into  F,  with  the 

k 

operator  norm 

IAI  ”  sup  {lA(v1#  ....  vk)lF»  Iv^lg  <  1,  i  -  1,  ...»  k> . 

I 

If  *:E  ♦  F  is  a  smooth  mapping,  its  k-th  Frechet  derivative  at  a  point  u  e  E  is 

|#  le 

D  4 (u)  e  L (E,F).  It  is  well  known  that  high-order  derivatives  are  symmetric 
v 

multilinear  mappings.  D  *(u)  is  therefore  completely  determined  by  assigning  its  values 

on  elements  of  the  form  v^«  (v,v,...,v)  e®  E.  We  write  B  for  the  closed  ball 

k  p 

centered  at  the  origin  with  radius  p,  and  T  ♦  for  the  n-th  order  Taylor  expansion  of 
the  map  4>  at  the  origin,  i.e. 

(T"*)(u)  -  \n  -rj  dVoJ’U111. 

i-0 

Given  a  function  *  »  ♦  (u,x)  defined  on  a  product  space  ExF,  its  partial  derivatives 
are  denoted  by  I"  this  case,  stands  for  the  n-th  order  Taylor  expansion  of 

4  at  (0,0).  The  Landau  order  symbols  O  and  o  will  also  be  used.  For  the  basic 
properties  of  differential  calculus  in  Banach  spaces,  our  general  reference  is  Dleudonne 
(2).  The  approximation  procedure  considered  in  this  paper  is  based  on  the  following 


simple  result: 


Proposition  1.  Let  e,f  be  Banach  spaces,  let  (u,x>  *  4(u,x)  be  a  Ck  tup  frosi  a 
neighborhood  of  the  origin  in  ExF  into  F  such  that  4(0,0)  «  0,  3^4(0, 0)  *>  0,  and 

let  the  map  u  :  E  ♦  F  be  implicitly  defined  by  4(u)  -  4(u,4(u)).  If  4,4  are  C* 
maps  such  that  Tn4  »  Tn4  and  Tn  \  ■  Tn  %  for  some  n,  0  <  n  <  k,  then  the  n-th 
order  Taylor  expansion  at  the  origin  of  the  map  u  ♦  T(u,4(u))  coincides  with  T°4 • 
Proof.  Since  the  maps  ♦  (•)  and  4(*,4(*)>  are  both  C*.  we  only  need  to  show  that, 
under  the  above  hypothesis, 

lim  I4(u)  -  4(u,4(u))l»lul  "  ■  0.  (2.1) 

n-»0 

By  Taylor's  theorem,  for  a  suitable  constant  C  and  for  u,x  sufficiently  close  to  the 
origin  one  has 

■4(u)l  <  Clul,  l*(u)l  <  Clul,  I4(u)  -  4(u)l  <  Clul", 

I3^4(u,x)l  <  Ct lul+lxl ) ,  l*(u»x>  -  4<u,x)l  <  C  (lul"+1  +  lxln+1). 

These  inequalities  imply 

■4(u)  -  4(u.4(u))l  <  I4(u,4(u) )  -  ft u,4(u))l  +  IT(u,4(u) )  -  4(u,4(u))l 

<  C(lu!n+1-H4(u)i”+1>  +  /’  l3jtf(u,54(u)+(  1-C >4(u> >1  *l4(u)  -  4(u>i<* 

<  C(1+Cn+1)luln+1  +  Cdul  +  Clul)  •  Clul”, 
which  proves  (2.1). 

Corollary  1.  If  4,  4  satisfy  the  assumptions  made  in  Proposition  1,  then 

lim  I  4(u,(Tn4)(u) )  -  (Tn4)(u>  Idul*"  -  0.  (2.2) 

u*0 

$3  The  Taylor  Formula 

Most  of  this  paper  is  concerned  with  the  system  (1.2).  Notice  that,  by  setting  uf  5  1, 
a  control  system  of  the  form  (1.1)  is  recovered  as  a  special  case  of  (1.2).  To  simplify 
all  further  discussion,  we  assume  that  the  vector  fields  g^  are  C  •  In  the 
following,  |g^(x)|  denotes  the  euclidean  norm  of  the  vector 


g^(x>  e  Rd,  while  l^g^(x)|  1b  the  operator  norm  of  the  dxd  matrix  of  first  order 
partial  derivatives  of  at  x.  The  set  of  admissible  controls  is  defined  as 

U  »  {u  “  (u.,  ....  u  )  e  (.'  (  [0,“)  >8?") )  |u.(t)|  <  1,  1  •  1,  ...,  m,  t  >  0). 

1  in  l 

The  L1  norm  on  the  set  of  controls  and  the  C°  norm  on  the  set  of  trajectories  will  be 
always  used.  We  write  (1.2)  in  the  more  compact  form 
x(t)  »  G(x(t))u(t)  ,  x(o)  “  0.  (3.1) 

G  is  therefore  a  C  map  from  Rd  into  L(rfn'R<5).  Call  G^(s)  its  j-th  derivative 
at  z  e  R?  and  let  GR  -  T"G.  Set  E  -  L1  ( (0 ,» )  i  R®) ,  F  -  C<  10,“ )  »Rd  > ,  and  define  the 
map  *  s  ExF  ♦  F  by 

•<u,y)(t)  ”  J*  G(y(s))u(s)  ds  .  (3.2) 

Then  •  «  o  with  ♦,(u,y)(t)  “  (u(t),  G(y(t))) 
and 

*"(u,Z)(t)  “  Z(s)u(s)ds. 

*'  is  a  C  substitution  operator  and  ♦“  is  bilinear.  Hence  *  la  k  times 
continuously  Frechet  differentiable,  for  all  k.  In  particular,  3^j  *(uo,xQ)(t)  is  the 
j-linear  map 

yfJJ  ♦  J*  G(^(x0(s))  yfJ,(s)  u0(s)  ds,  (3.3) 

3  3^  *(u  ,x  )(t)  is  the  multilinear  map 
U  X  o  o 

(u.y111)  ♦  J*  G(i,(x0(s))  y{J)(s)  u(s)  ds  (3.4) 

and  3*3^  *  5  0  for  i  >  1,  because  the  dependence  on  u  is  linear.  Notice  that  the 
input-output  map  u(*)  *  x(u,*)  «  $<u)(»)  generated  by  (1.2),  or  equivalently  by  (3.1), 
is  implicitly  defined  by  the  equation  f(u)  »  $(u,+(u)),  and  that  both  *  and  3x  * 
vanish  at  (0,0).  The  Taylor  expansion  of  ♦  at  the  origin  can  therefore  by  computed 
recursively,  by  means  of  Proposition  1.  We  will  show  that  Tn*  can  always  be  written  as 
a  finite  sum  of  certain  iterated  integrals  of  the  control  u,  here  called  integral 
monomials  (integromials,  in  short). 

Definition.  The  family  M  ■  M(m,d)  of  integral  monomials  is  the  smallest  set  of  mappings 
v  from  i  ( (0,«)  i r"*)  into  Rd  with  the  following  properties: 


-4- 


(3.5) 


i)  For  every  linear  A  »  rf1  ♦  Rd,  the  map 
U  «  (u,t)  ♦  J*  A«u(e)ds 
ie  in  M 

il)  If  |if>  ...»  |tk  e  M  and  if  B«(Rd)k  ♦  ie  k-linear,  then  M  also  containa 

the  sap 

V  I  (u,t)  ♦  J*  B(u,(u,s),  ....  U^tu.e))  u(e)de.  (3.6) 

Using  a  canonical  identification,  we  shall  regard  integral  monomials  as  multilinear 

mappings  frow  L1((0,»)iKB)  into  C°( [0«**)»Rd).  If  V  6  M  is  k-linear,  we  say  that  W 

has  order  k.  Notice  that  if  in  (3.6)  has  order  v^i  -  1,  ...,  k),  then  |i  has 

order  1  +  +  ...  +  vfc.  Consider  now  the  approximate  system 

x(t>  -  Gv(x(t))u(t)  ,  x(0)  -  0,  (3.7) 

recalling  our  definition  i  Gv  ■  TUG.  the  first  Picard  iterate  for  (3.2) 

P,(u,t)  “  ffc  G(o)u(a)ds  (3.8) 

1  o 

is  an  integral  monomial  of  the  form  (3.5).  Moreover,  if  the  n-th  Picard  iterate  PR  can 

be  written  as  a  finite  sume  of  integral  monomials,  say 
N(n) 

P  (u,t)  -  £  |i,(u,t), 

1-1  J 

then  the  same  holds  for  Pn+1.  Indeed 

P  ,(u,t)  -  J*  IV  G<k)(0)(P  (u,s)),K1u(s)ds 
>  k-0 

can  be  reexpanded  into  a  finite  sum  of  integral  monomials  of  the  form  (3.6),  namely 

Pn+1(u't)  ■  r*  r  ll  jJ-  G(k)<°)  . Wo{k)(u,s))uts)ds,  (3.9) 

km0  <ter 

where  T  »  T(k,N(n))  is  the  set  of  all  mappings  a  from  {1,2, ...,k)  into 
{  1, • « • ,N(n)} . 

In  practice,  the  k-th  order  Taylor  expansion  T*+  at  the  origin  for  the  input-output 
map  +!u(<)  *  x(u,*)  generated  by  (3.1)  can  be  obtained  by  either  of  the  following 
methods  s 

-5- 


I)  Compute  the  k-th  Picard  iterate  for  the  system 

x(t)  -  Gk.1  (x(t))u(t)  ,  x(o)  -  0  (3.10) 

and  discard  the  integral  monomials  of  order  >  k. 

II)  Set  *n(a,,l  “  (Tnip  )(u)(«)  and  recursively  derive  xn  from  (n  “  1,...,k)  by 

expanding  the  map 

<u,t)  ♦  J*  Gn_1(xn_1(u#s))  u(s)ds  (3.11) 

into  a  sum  of  integral  monomials,  discarding  those  which  have  order  >  n# 


Indeed,  for  n  >  1  set 

*n<u,x)(t)  «  Gn_t(x(s))  u(s)ds  (3.12) 

and  check  that  the  Taylor  expansions  at  (0,0)  of  *n  and  *  (defined  at  (3.2))  coincide 


up  to  order  n.  Let  Pn  be  the  n-th  Picard  iterate  for  (3.10),  regarded  as  a  map 
u(*)  +  Pn<u,  • )  from  L1  into  C°.  Assume  that  Tn  1p  .  !MTn~1<|>.  Then,  by  setting 


f  ”  *  “  pn_i»  Proposition  1  yields  Tnpn  “  T*V ,  provided  n  <  k.  By  induction. 

Observing  that  pk  is  a  finite  sum  of  integromials,  to  obtain  its  k-th 


order  Taylor  expansion  at  the  origin  one  merely  discards  the  terms  of  order  >  k.  This 
justifies  I).  Now  set  T  ”  *n,  <|>  ”  Tn  %•  By  Proposition  1  the  map 

«  *  ^(u^T"  %)(u)),  otherwise  defined  at  (3.11),  has  the  same  n-th  order  Taylor 


expansion  at  the  origin  as  <f.  Discarding  from  (3.11)  the  integral  monomials  of 
order  >  n  we  thus  recover  from  $  ■  One  obtains  by  repeking  the  above 


procedure  for  n  «  1,  ....  k. 


§4  Examples 

Example  1.  The  third  order  Taylor  expansion  for  the  scalar  system 
x ( t )  *  cos  x(t)u1(t)  +  u2<t),  x(o)  -  0 

is 

<T3*Hu)<t)  -  ft(u,(s)  ♦  u ,  (s)lds  -  V7  u  (0,  )  (  /  1u,(o,)dc.)2dc, 

Oi  2  0110122  1 
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If  uj  is  constrained  to  be  identically  1  we  obtain  an  approximation  for  the  system 
x(t)  *  cos  x(t)  +  u( t ) ,  x( 0 )  -  0,  namely 
x(u,t>  -  t  -  +  VfeJ*  u(<J)do  -  V2/£  (  /°1u(o2)d<J2)2d®1 

t  °1  4 

~  Ja  Jo  u(o2)do2do1  +  0(t  ). 

Example  2,  Consider  on  R3  the  control  system 
x  »  X(x)  +  Y(x)*u  ,  x(0)  -  0.  (4.1) 

Assume  that  X(o)  -  0  and  span  f(advX,Y)i  v  »  0,  1,  2)  -  R3.  then  the  map  *  «  <sv 

S2,  S3)  ♦  (exp  S1  Y)  •  (exp  s2(y,X)>  .  (exp  s3r[Y,X],X])  (o)  defines  a  local  chart  of 

a  neighborhood  of  the  origin  in  R3.  In  this  chart,  the  third  order  Taylor  approximation 
for  (4.1)  takes  the  form 


x^u.t)  -  /*u(s)ds  +  /*  (  ]g 1  u(<»2)dc2>2do1  +  0(t4) 


k  o 

x2(u,t)  -  I*  jo1u(o2)do2doi  +  ~  J*  (  /O1u(o2)do2)2do1  ♦  0(t4) 


x3(u.t)  -  J‘  jj  /o2u(03)d<J3  dOjdo,  ♦  -§  /*  (  j\(t,2)ao2)2dai  ♦  0(t4), 

where  the  constants  kj,  k2,  k3  are  obtained  from  the  linear  relation 
k,Y(o)  +k2[Y,X](o)  +  k3[(Y,X],X](o)  -  [Y,[Y,X]J<o). 

To  prove  (4.2),  we  need  to  compute  some  Taylor  coefficients  of  X  and  Y  in  the  given 
chart.  The  definition  of  x  Implies  that  Y  =  (1,0,0),  [Y,X](o)  -  (0,1,0)  and 
t[T«X],  X] (o)  "  (0,0,1).  Re  thus  have 


~  (o)  -  Y«X(o)  -  (Y.Xl(o)  -  (0,1,0) 
3“l 


(o)  -  IY,X]  •  X(o>  -  l(Y,X],X](o)  -  (0,0,1) 


1  and  use  Theorem 


(O)  -  Y(Y»X) (o)  -  [Y,[Y,X]1(0)  -  (k1,k2,k3). 

3S1 

One  can  now  multiply  the  vector  field  X  by  a  fictitious  control  uq  5 
1  to  get  (4.2). 

Remark:  It  is  clear  that  the  truncated  Taylor  expansions  are  not  invariant  uni*  change 
of  coordinates.  The  above  example  shows  how  the  local  Lie  algebraic  structure  the 
system  can  yield  a  chart  in  which  the  expansion  takes  the  simplest  possible  fr  For  a 

general  procedure  to  construct  "canonical"  charts,  see  [S  p.  127]. 

is  Error  bounds. 

High-order  Taylor  expansions  are  the  primary  tool  in  the  local  study  of  functionals  at 
points  of  singularity.  If  the  origin  is  a  stationary  point  for  a  map  A  :  E  ♦  R,  one 

classically  proves  that  A  attains  a  local  minimum  at  0  by  showing  that  for  some  n  >  1 

(TnA  )  (u)  -  A  (o )  »  C  •lul"  (5.1) 

for  all  u  is  a  neighborhood  of  the  origin.  One  would  hope  to  use  a  similiar  argument 
and  prove  non-controllahillty  results  concerning  the  system  (1.1).  For  example,  let  f(o) 
•  0.  Then  (1.1)  is  not  locally  controllable  at  the  origin  if,  for  some  nontrivial  z*  @ 

T  >  0,  the  functional  A(u)  •  <  t*,  x(u,T>>  attains  a  local  minimum  at  u  S  0. 

Coercivity  conditions  such  as  (5.1),  however,  can  seldom  be  obtained  in  connection  with 
control  systems.  To  overcome  this  difficulty,  one  needs  bounds  on  the  error 
c(u)  ■  l(T%)(u)  -  4(u)l  which  are  sharper  then  the  classical  bound  e(u)  <  Cluln+1 
supplied  by  Taylor's  theorem. 

Lemma  1.  Let  x(u,*)  be  the  solution  of  (1.2)  corresponding  to  the  control  u.  Then 
there  exists  a  tQ  >  0  such  that,  for  every  continuous  map  y  from  [0,tol  into  the  unit 
ball  B,  £  ifi, 
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(5.2) 


I 


|y(t>  -  x(u,t)j  <  2  sup  (I  /*  G(y(s) )u(s)ds  -  y(x)|»  0  <  T  <  t} 
for  all  t  e  [0,to],  u  e  0. 

Proof.  Sat  t0  “  M_1 ,  with 

N  .  2b  •  SUp  {jgi(x)|  +  |7g1(x)|,  x  e  B^,  1  -  1,  ....  m)  . 

Than  the  nap  *  dsflned  at  (3.2)  satisfias 

'  ♦(u'yi>  -  ♦(u'y2),Cto,t  )  <1/2'yi  ‘  Vcio.t  ] 

o  o 

and  *(u,y1>(t)  t  1^,  for  avary  u  e  U,  t  e  [0,to]  and  every  continuous  maps 
y1«  y2  :  tO.tl  ♦  .  Since  x(u,*)  is  a  fixad  point  for  tha  may  y  *  *(u,y),  the 

contraction  napping  theorem  yields  I y  -  x(u,*)l  <  2ly  -  d(u,y)l,  i.e.  ,  (5.2). 

tasau  2.  For  each  u  e  v,  let  x(u,*)  be  the  solution  of  (3.1)  and  let  yn(u,» )  be  the 
solution  of 

y(t)  -  Gn_t(y(t))u(t),  x ( o )  -  0.  (5.3) 

Than  there  exist  C  and  t  >  o  such  that 
|yn(u,t)  -  x(u,t)|  <  C  |  J*  r.(o)u(a)do!Bds  (5.4) 

for  all  t  e  (0.to] ,  u  e  0. 

Proof.  Choose  t,  >  0  so  small  that  yn(u,t)  e  B1  for  all  u  e  U,  t  e[0,t.,].  Lemma  1 

implies  the  existence  of  a  t  >  0  such  that  t  <  t.  and 

|y  (u,t)  -  x(u/t) |  <  2  sup  {  |  j*  G(y  (u,s))u(s)ds  -  y  (u,T)|,  0  <  x  <  t)  (5.5) 

n  o  n  n 

for  all  t  €  (0,toJ ,  u  «  0.  Using  (5.3)  we  have 

|y  (u,t)  -  x(u,t)|  <  2  sup  {  |  / T [G(y  (u,s))  -  G  (y  (u,s) )) u(s)da I ,  0  <  x  <  t} 
n  on  n- 1  n 

<  2  •  C1  |yn(u,s)|nds,  (5.6) 

where  tha  constant  C,  depends  only  on  the  size  of  the  n-th  derivative  of  G  on  B1 . 

Set 

M  »  n  •  sup  {|V(Tn1gi)(z)|izeB1^m<1,i“1,...,m>,  | 

n(t)  »  |  j*G(o)u(s)ds|  ,  w(t)  -  |yn(u,t)  -  j^G( o)u( s ) da |  .  For  almost  every  j 

t  «  [0,t  J,  w(t)  <  I  G  ,(y  (u,t) )  -  G  (o)  !  <  M  •  |y  (u,t)|  <  M(w(t)  +  n(t>).  I 

o  n* »  n  n-i  n  | 


t  !TC_ 

Since  w(o)  *  0,  Gr on wall’s  lemma  yields  v(t>  <  ]Q  n(s)ds#  with  C2  ■  M*  •  Using 

M 

Holder's  inequality  we  get  the  estimate 


2  /*  |y(u,s>  |nds  <  2  /*  fn(T)  +  C2  n(a)ds) ndt 

<  J*  (2ri(T))ndt  +  J*  (2C2  J*  n(s)de)ndT  (5.7) 

<  jt  (2n(T))ndT  +  ( 2c  )n*t  (  J*  nn(s)ds)tn_1  <  c.  jt  nn(r>dr, 

O  /  O  J  O 


with  C3  «  2n  +  (2C2  tc)n.  Combining  (5.6)  with  (5.7)  one  has  (5.3),  with  C  ■  C^Cj. 

The  above  lemmas  provide  a  simple  estimate  for  the  error  in  the  Taylor  expansion  of  $ . 

Theorem  1.  Let  x(u,*)  he  the  solution  of  (3.1).  Set  Gn_i  “  Tn-1G, 

x  (u,t)  •  (Tn$ ) ( u) ( t) .  Then  there  exist  C,  t_  >  0  such  that 

n  0 


|x(u,t)  -  XR(u,t)  !  <  C  I  j®  G(o)u<0>  do|nds  + 

+  2  sup  {  |x  (u,t)  -  JT  G  ,(x  (u,s) )u( 8>ds ] »  0  5  T  5  t} , 
n  o  n- 1  n 


(S.8) 


for  every  u  e  U,  t  e  [0,to) . 

Indeed,  |x(u,t)  -  x^U/t))  <  |x(u,t>  -  yn(u,t)|  +  |yn(u,t)  -  x^u.t)].  The 

estimates  (5.4)  and  (5.2)  with  G  replaced  by  yield  (5.8). 

d  2  3 

Example  3.  Consider  the  bidimensional  system  —  (x,y)  “  (u+y*cos  x,  x  -x  ), 
(x(0),y(o))  «  (0,0).  Its  third  order  approximation  is 
x(u,t)  *>  u(s)d8  +  e^lu.t)  ,  y(u,t)  «  (  /*  u(o)do)2ds  +  e2(u,t)  . 

Theorem  2  yields  the  bounds 

o  o 

|ei  (u,t)  |  <  C  1  /*  u(o)do|3ds  +  2  J*  jo  <  /o2u(o3)do3)2do2do1,  (5.9) 
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for  i  »  1,  2,  u  e  U  and  t  sufficiently  snail*  From  (5.9)  we  see  that 
(c2(u,t)|  *  o  (  /*  (  /*  u(<J)d<r)2ds),  hence,  for  snail  t,  x(u,t)  >  0  and  the  ays  ten  is 
not  locally  controllable.  For  any  fixed  tQ  >  0,  consider  the  control  (t)  ■  cos  It. 

As  X  ♦  •,  lul*  tends  to  (2  t  A)*»  while  Jt  (  /*  u,  (o)(do)2ds  tends  to 

iVt  ]  °  o  °  x 

zero.  Bounds  of  the  type  I  e^(u,t)|  <  Ct  or  |e^(u,t)|  <  C  • I ul  are  therefore  too 
weak  for  proving  non-controllability  even  in  this  sinple  case. 

{6  Directional  Estimates. 

To  obtain  more  precise  bounds  on  the  error  in  the  Taylor  expansion  of  (1.2),  in  this 
section  we  split  tP  into  a  sum  of  orthogonal  subspaces  Vp  and  eatinate  the  else  of  the 
error  separately  on  each  Vp  Given  the  control  system  (1.2),  define  an  increasing 
sequence  of  subspaces  Sp  R?5  recursively  by  setting  i)  80  “  (0),  ii)  Sp  is  the 

smallest  subspace  of  if1  such  that  for  all  i  “  1,  ...,  n  and  k  “  0,  . ..,  p  one  has 

Dk  gA(o)  (z1# . . .  ,*y )  e  Sp  for  every  k-tuple  (c^,...,*^)  with  s^  e  and 

♦  J2  +  . . .  +  jk  <  p.  In  particular,  S ,  is  the  smallest  subspace  that  contains  the 

m  vectors  g1(o)  and  is  Invariant  under  the  linear  operators  Vg^o)  (i  “  1,  ...,  m). 

Now  choose  any  p  >  1.  For  1  <  p  <  p  define  Vp  as  the  orthogonal  complement  of  Sp_j 

in  S  .  Let  V_  be  the  orthogonal  complement  of  8  in  Finally,  denote  by  x 

P  d  P-'  P 

the  orthogonal  projection  R  ♦  .  With  this  notation,  we  have 

Theorem  2.  Let  u(*)+x(u,*)  be  the  input-output  map  generated  by  (1.2)  and  let 

u(.)  ♦  xn(u,.)  be  its  n-th  order  Taylor  approximation  about  the  null  control.  If 

p  e  (l,...,p)  and  n  >  p,  then  there  exist  CQ,  tQ  >  0  such  that 

I x  (x(u,t)  -  X  (u,t))|  <  C  tn_p  /*G(o)u(d )do |Pda  (6.1) 

p  n  o  o  o 

for  all  t  e  (0,to),  u  e  V. 

Proof .  Let  yn  be  the  solution  of  (5.3).  By  Lemsui  2,  the  difference  yn(u,t)  -  x(u,t) 
satistifies  a  bound  of  the  form  (6.1)  for  any  p  -  1,  ...,  p.  Therefore  we  only  need  to 
show  that 
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where 


|«  (y  (u,t)  -  x  (u,t))|  <  C  tn  P  /  |x  (u,s) |pds, 
p  n  n  o  1 


x^u.e)  «  /  G(o)u(O)d0 
o 

gives  the  first  order  Taylor  approximation  to  the  trajectory  x{*,»).  In  §3,  xn  was 
proven  to  be  a  finite  sum  of  integral  monomials  occuring  within  the  first  n  Picard 
iterates  for  the  system  (5.3).  To  keep  track  of  its  size,  to  each  one  of  the  above 
integromials  we  attach  an  integer  Y(li)  as  follows:  If  u  -  x^,  defined  at  (6.3), 
set  Y (u )  “  1.  If 

u(u,t)  -  J*  G(lc>(o)(u1(u,B),...,uk{u,s))u(8)da,  (6.4) 

set  Y(U)  “  Y(U1)  +  •••  +  Y(Uk>.  We  stress  that  this  definition  refers  exclusively  to 
the  integromials  arising  via  Picard  iterations  for  the  particular  system  (5.3)  presently 
considered.  Comparing  the  definitions  of  Y  and  of  the  subspaces  Sp,  it  is  clear  that 
U(u,t)  e  for  *11  t  >  0,  u  e  0.  Also  notice  that  if  the  order  of 

U  is  v  >  1,  then  Y(U)<v-1.  A  basic  estimate  on  the  Bize  of  integral  monomials  is 
now  given. 

Lemma  3.  Let  u  be  an  integral  monomial  of  order  v  >  1  occurring  in  some  Picard 
iterate  for  (5.3),  and  let  Y(u)  “  p.  Then  there  exists  a  constant  C  such  that 
I U (u ,t ) |  <  C  tV  P~1  |x1(u,s)|Pds  (6.5) 

for  all  u  e  U,  0  <  t  <  1. 

Proof.  Notice  first  that  for  the  integromial  x^*,*)  ire  have  v  “  p  -  1.  For  all 

others,  v  >  1.  Moreover,  the  only  integral  monomials  for  which  v  -  p  +  1  have  the 
form 


p(u,t>  «  J*  ■V  C^No)  (x  (u,s))^u(s)  ds  (6.6) 

O  pi  I 

and  clearly  satisfy  an  estimate  of  the  type  (6.5).  The  general  case  will  be  proved  by 
induction  on  v,  assuming  v  >  p  +  2.  If  u  is  given  by  (6.4),  let  have  order 

<  v  and  let  y(u^)  -  pt  (i  *  1,...,k).  The  inductive  hypothesis  implies  that  either 


|Ui(u,s)l  <  C  • 


Vpi"’ 


/"  I x f  ( u ,o )  t  1  do 


(6.7) 
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(or  mm  conatant  and  all  u  €  U,  0  <  t  <  1. 

Ualnq  tha  ayaaetry  of  G*k *(o),  ve  can  raordar  tha  u 1  and  aasuaa  that  (6.4)  haa  tha  fora 

g(u,t)  -  J*  jr  G<lt,(o)(u1(u,a),...,uh(u,a),  x, (u.a) , . . .  ,x1  <u,e)  )u(a)da  (6.8) 

(or  aoat  0  <  h  <  k,  with  aach  g^(1  <  i  <  h)  having  ordar  >  t.  Tha  ordar  of  g  in 
(6.8)  la  than  v  ■  1  ♦  v1  ♦  •••  ♦  vh  ♦  k  -  h,  while  p  "  T(g)  ■  p,  ♦  •••  p^  ♦  k  -  h.  In 

tha  following,  va  aat  q  -  p,  ♦  •••  ♦  ph  -  p  -  k  ♦  h.  If  q  -  0,  than  w  haa  tha  for» 

(6.6)  and  tha  estimate  (6.5)  ia  iaaadiate.  If  q  >  0,  tha  inductiva  hypothaaia  and 


■ 

Noldar'a  inaquality  imply 

|g(u,t) I  <  c'  J*  |x  (u,a)!k'h  n  (c  a  1  Pi  J* 

i-1 


<  C  j*  . ‘v-u-p ^  [(]•  i*1(«.B)i'«d»Pl/ 


(q-p, )/q 

a  )  da 


<  c  |t  av-p-2  (x  (tt,a)|k-h(  /•  |x,(u.o)|qdB)  da,  (6.8) 

O  '  o  » 

whara  tha  oonatanta  C,  C  ,  ara  indapandant  of  u  and  t.  If  h  -  k,  (6.5)  ia  a 
trivial  conaaquanca  of  (6.9).  If  h  <  k,  ona  racovara  again  (6.5)  froai  (6.9)  integrating 
by  parta  and  uaing  Holdar'a  Inaquality i 

(g(u.t)l  <  c  tw‘p'2(  /*  (x,(u.a)|k'hda){/tfx,(a,a)l<,da) 

O  1  O  1 

<  c  tw-P'1  St  |x.(u,a)|k“h'*qda. 

O  * 

Returning  to  tha  proof  of  Diaoraa  2,  tor  any  control  u  e  O  and  any  conatanta 

C,  t  >  0,  datlna  ),tt<C,t)  aa  tha  aat  of  all  nape  a  e  C([0,»)i  Rd)  auch  that 

|w  (a(t)  -  x  (u, t ) ) I  «  c  •  tn_p  /*  (x  (u,a)|pda  (6.10) 

p  n  o  * 

tor  all  t  e  [0,t],  p  •  1 . .  We  data  that  y  (u,«)  e  )  (C  ,t  )  for  aome  C„, 

n  u  o  o  o 

tQ  >  o  and  all  u  e  0. 

Laaaa  4.  there  axiat  conatanta  CQ,  tQ  auch  that,  for  avary  u  e  0,  tha  nap  9u  dafinad 
by 

9  ( x ) ( t )  -  Jfc  G  , (a(a) )u(a)da 
u  o  n- 1 

mapa  ) (C  ,t  )  into  itaalf. 
u  o  o 

a  - 

Proof.  Given  u  e  0,  -  '  C((o,«)i  R  ),  for  p  *  1,  ...»  p  aat 
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w  (t)  -  *  (r(t)  -  x  (u,t)) 
P  P  n 


We  than  hava 


l»  (x  )(u,t)  -  *  <*<t)  |  <  |w  (x  (u,t)  -  Jt  G 
p  n  u  p  n  4  o  n- 

♦  !*  (  j*  G  (x  (u,a))u(a)da  -  Jt  G  (x  (u,e) 
p  Jo  n-1  n  Jo  n-1  n 


(x  (u,e) 
i  n 


>u(s)da)  I 

(a) )u(s)da) I 


-  I A(u,t ) I  +  |B(u,t)|.  (6.11) 

Notice  that  A(u,t)  is  a  finite  sun  of  integral  monomials  which,  by  Corollary  1  in  {2, 

have  order  >  n.  If  u  is  one  of  such  monomials,  then  either  y(u)  >  p  or 

*  (u(u,t) )  I  0.  By  Lemma  3,  there  exists  a  constant  C.  >  0  such  that 

P  1 

I  *  (A(u,t) ) I  <  C,  tn'p  jt|x,(u,s)|pde  (6.12) 

P  1  O  I 

for  all  u  e  U,  0  <  t  <  1,  p  »  1,  ....  p".  Take  CQ  -  2C1.  To  determine  a  suitable  tQ, 
observe  that  B(u,t)  can  be  written  as  a  finite  sum  of  terms  of  the  type 


A(u,t)  -  J*  G(k)(0)  (wp  (s),  ....  wp  (s),  ut+1(u,s),  ...,  yl+h(u,s), 
X^u.s),  ...,  X  ^  (u  ,  B  )  )u(  S  )  ds  , 


t  >  1  and  the  integromials  ut+1,  ...»  hava  order  >  1. 

+  h,  let  p^  ■  yty^).  Observe  that  the  order  of  is  then  at 

For  i  »  t  +  1 

. 1 

least  pt  ♦  1. 

if  1  e 

),u(Co,  t)  we  thus  have  the  inequalities 

Iut(u,s)l  <  C 

1“  ,x1(u' 

P< 

,a)|  do, 

|w  (u, s)  |  <  C 

.""P1  J 

f"|x  (u,o)|  1do, 

O  1 

Pi 

0  ' 

for  a  suitable  constant  C  and  s  <  min  {t,1}.  The  definition  of  8^  implies  that 
either  p^  ♦  ...  ♦  P^+h  +k-i-h>p,  or  »p(A(u,t))  =  0.  The  same  arguments  used  in 
the  proof  of  Lemma  3  now  yield  constants  C2*  tj  >  0  such  that 
*  (A(u,t) )  <  C,  tn_p+1  /t|x,(u,s)|pds,  (6.13) 

P  2  Ol 
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<6.14) 


for  0  <  t  <  tj,  u  e  U. 

Therefore,  there  axiat  C3,  t3  >  0  such  that 
|»  (B(U't) ) I  <  C,  tn_p+1  /‘lx.tu.sJlPds 
for  all  u  e  0,  p  ■  1,  ....  p,  t  8  [0,t j] .  Comparing  (6.11)  with  (6.12)  and  (6.14),  wa 
saa  that  Leans  4  holds  with  tQ  -  min  {1,  tj,  C^Cj'). 

The  conclusion  of  the  proof  of  Thaorea  2  is  now  straightforward.  For  all  u  6  0, 
yn<u,*)  is  the  unique  fixed  point  of  0U«  By  Leona  4,  yn(u,*)  6  (  cQ « tQ ) ,  hence 
(6.2)  follows  from  (6.10). 

With  the  saae  notation  of  Theorem  2  we  have 

Corollary  3.  If  p  >  1,  n  <  p,  then  thera  exist  CQ,  tQ  >  0  such  that 

I  »  (x(u,t)  -  x  (u,t))|  -  |  *  (x(u,t))|  <  C  /*|  /"  G(o>u(C)dc|pds  (6.15) 

p  n  p  o  o  o 

for  all  u  e  0,  t  e  [0,tQl,  p  “  1,  ...»  p. 

Indeed,  x,,  is  a  sua  of  integral  aonoaUala  having  order  <  n.  Hence 

Y(ut)  <  p  and  vp(wi(u,t))  =  0.  This  iaplies  Xplx^u.t))  s  0.  Setting  n  «  p, 
Thaorea  2  yields  the  bound  (6. IS). 

Kxanqjle  4.  consider  on  II3  the  system  x  -  <xf,  x2,  x3>  -  (u  cos  x,  -  x2  -  x3, 

,23  , 

tg(x,  -  x3),  sin*  x,  -  x2>  x,),  x(0)  ■  0  «  r.  »  third  order  expansion  yields 

.  * 

x,(u,t)  -  r  u(s)ds  -II  I  U(0,)d9  do,*,  *  C.(u.t), 

1  O  o  O  O  j  J  4  I  1 

..  ff. 

Xj(u.t)  -  J*  Jo  u(o2)do2do1  ♦  c2(u,t)  , 
x  (u,t)  *  ]t  (  /*  u(o)do)2ds  ♦  c,(u,t). 

3  O  O  J 

Tor  this  systea,  81  -  V1  » 

{(C,,  C2,0)  ;  *  **'  s2  “  *3'  v2  “  ^°»0»53)»  t3  «  *)•  By  Thaorea  2  there  exist 

C,  T  >  0  such  that 
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